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Abstract. By presenting linear and nonlinear gyrokinetic studies, based on a 
balanced neutral beam injection deuterium discharge from the DIII-D tokamak, we 
demonstrate that impurities alter the scaling of the transport on the charge and mass 
of the main species, and even more importantly, they can dramatically change the 
energy transport even in relatively small quantities. A poloidally varying equilibrium 
electrostatic potential can lead to a strong reduction or sign change of the impurity 
peaking factor due to the combined effect of the in-out impurity density asymmetry and 
the E x B drift of impurities. We present an approximate expression for the impurity 
peaking factor and demonstrate that impurity peaking is not significantly affected by 
impurity self-collisions. 



1. Introduction 

Several models have been proposed to explain the favorable isotope mass scaling of 
the energy confinement in tokamaks [H EJ [3j HJ [5], yet a commonly accepted and 
robust theoretical explanation of the isotope effect which is consistent with the other 
experimentally established parameter scalings [HI El El E] is still lacking. When the 
turbulent energy transport is determined by the properties of the main ion species, 
dimensional considerations suggest that the charge and mass scaling of the heat 
diffusivities, x, should follow the gyro-Bohm scaling, XgB ~ Pi v i/ a ^ y/^/Zf, where 
Vi = (2T i /mj) 1 / 2 is the ion thermal speed, a is the plasma minor radius, pi = Vi/Vt ci 
is the ion thermal Larmor radius, with the cyclotron frequency Q C i = ZieB/rrii of an 
ion of charge number Zi and mass nij, and B is the magnetic field strength. Studies of 
charge and mass effects on the turbulent energy transport in DIII-D tokamak plasmas 
show that the presence of impurities can contribute to deviations from a pure gyro- 
Bohm scaling [TT] . In this paper we use the expression "pure gyro-Bohm" if differences 
in the turbulence between similar plasmas with different main species can be exactly 
transformed out by normalizing temporal and spatial scales to a/vi and pi, respectively. 
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In the first part of the paper we present similarity studies, that is, a series of gyrokinetic 
(GK) simulations with gyro [10 j where the local geometry and profile parameters - 
taken from a deuterium discharge studied in Ref. [TTJ - are kept fixed while the ion 
composition is artificially changed from case to case. First the effect of impurities (in 
particular, carbon or deuterium) on energy transport in hydrogen isotope and helium 
plasmas is studied. Then we investigate the role of the main species on impurity 
transport, and conclude that impurity transport does not follow the pure gyro-Bohm 
scaling (based on the properties of the main ion species) so closely as the energy 
transport. 

One of the ways to reduce impurity accumulation in tokamak cores is to apply 
central radio frequency (RF) heating [TJl EH [EH US] . RF-heating can generate poloidal 
impurity density asymmetries [16] and if the asymmetry is sufficiently strong, it can lead 
to a sign change in the impurity peaking factor [T7"t IT8] . Recently it was shown that a 
sign change of the impurity peaking factor can happen at much weaker (that is, realistic) 
asymmetry strengths than were considered in Ref. [18] if the E# x drift of impurities 
in the poloidally varying equilibrium electrostatic field is taken into account [TS]. In the 
second part of this paper, we extend the work of Ref. [19], by including parallel dynamics 
and finite Larmor radius (FLR) effects and modeling impurity self-collisions with the 
full linearized operator for impurity self-collisions to derive an improved approximate 
expression for the impurity peaking factor. 

2. Charge and mass effects on turbulent transport 

In this section we study how the charge and mass scalings of the energy transport 
are affected by the presence of impurities, through linear and nonlinear GK similarity 
studies with gyro, based on local profile and magnetic geometry data from an L-mode 
phase of a weakly rotating DIII-D deuterium discharge (129135, 1250-1300 ms). For 
detailed information on the resolution of the simulations and on the studied discharge 
we refer the reader to Ref. [UJ. The dominant impurity in the discharges studied 
in Ref. [JTJ was carbon, and the hydrogen and helium discharges had a significant 
deuterium fraction present; this motivates our choice of impurities in the simulations 
presented in this section. We neglect particle flows within the flux surface and to model 
different ion compositions we keep the electron density profile fixed and assume constant 
ion concentrations. For different main species with the same charge the results of local 
GK simulations - normalized to main ion species units - should exactly coincide when 
impurities and collisions are neglected and adiabatic electron response is assumed [TT] . 
Thus we normalize growth rates, 7 to c s i/a, where c S i = (Te/nii) 1 ^ 2 , binormal wave 
numbers k y to 1/ p S i = Q C i/c S i, and fluxes to gyro-Bohm units Q 9 b% = n e T e c S i(p S i/a) 2 , 
FgBi = n e c si (p si /a) 2 oc m\ l2 /Z 2 . 

In Figs. [TJ (a,b) the growth rate of linear modes at k y p S i = 0.3 (at r/a = 0.55) 
is shown as a function of carbon concentration nc/n e . The local geometry and profile 
parameters, in the notation of Ref. [10], are: q = 1.78, s = 0.75, k = 1.34, 5 = 0.11, 
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a/L n = 0.68, a/L Ti = 1.89, a/L Te = 2.69, T;/T e = 0.82, v ei = 0.14 c s /a. The solid 
line corresponds to a deuterium plasma, the dashed is a hydrogen and the dotted 
is a helium plasma. We also show simulations for diluted deuterium plasma (long 
dashed lines), i.e., the deuterium concentration is decreased to the value that would 
correspond to the given carbon concentration (according to rii/n e = 1 — Zc<nc/n e ), 
but the carbon has no response to fluctuations (the carbon density is set to zero in 
the simulation). Fig. [TJ (a) shows the simulation results with adiabatic electrons. The 
decrease in the growth rates with increasing carbon concentration is practically linear. 
For lower carbon concentrations the slope of the 7 curves is the same in main ion 
species units, independently of the mass or charge of the main species. For higher 
carbon concentrations a small difference between hydrogen isotopes and helium appears 
regarding the slopes, but the deuterium and hydrogen curves still overlap almost exactly, 
in spite of the difference between the mc/rrti mass ratios. Noticeable deviations between 
growth rates in different hydrogen isotope plasmas appear only at higher wave numbers; 
noting that l/p s c ~ 2.45/p s ;. 

Figure H](b) shows trapped electron (TE) mode growth rates in a collisionless drift- 
kinetic electron simulation. The TE modes are weakly destabilized as the impurity 
concentration is increased and this effect is somewhat stronger in hydrogenic plasmas. 
This destabilization is due to the increased weight of trapped-electrons in the Poisson 
equation due to the diluting effect of impurities (similarly, the stabilizing effect of 
impurities on ITG modes mentioned before is mostly due to the decreased weight of 
ions). But it has to be noted that these TE modes are strongly stabilized when electron- 
ion collisions are introduced in the simulation, and the stabilizing effect of increasing 
Z e g appearing in the electron-ion collision frequency can be much stronger than the 
destabilizing effect observed in the collisionless case. 

Finally, we compare the deuterium plasma simulations with the correct physics 
(solid lines) to the model where only the diluting effect of impurities is taken into account 
(long dashed lines). The presence of carbon contributes to a further stabilization of the 
ITG modes, and this contribution is comparable to the effect of the dilution of the main 
species. However, for the TE mode, the presence of carbon has almost negligible effect 
compared to that caused by the dilution of the main species. 

Moving on to nonlinear simulations, the curves in Figs. Q] (c-f) represent the 
distribution of energy (or particle) fluxes over poloidal wave numbers; the kep S i integral 
of the curves give the fluxes in gyro-Bohm units; the corresponding flux values are 
given in the plot legends. The nonlinear simulations include kinetic electron response 
and electron-ion collisions and they are based on the local parameters of the studied 
deuterium plasma at r/a = 0.65 (q = 2.08, s = 1.1, k = 1.36, 5 = 0.14, a/L n = 0.79, 
a/L Ti = 1.94, a/L Te = 3.07, T.jT, = 0.91, u ei = 0.23 c s /a). 

First, we study the respective and combined effects of the carbon and deuterium 
impurities on energy transport. Figured] (c) shows the total ion energy flux in species 
units for a pure hydrogen plasma (solid curve), in hydrogen plasmas with carbon 
(dashed) and deuterium (dash-dotted) impurities, and with both carbon and deuterium 
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Figure 1. (a,b) Growth rates of the k y p S i = 0.3 mode as a function of carbon 
concentration. Solid: main ion species is deuterium, dashed: hydrogen, dotted: helium, 
long dashed: diluted deuterium (no carbon, non-quasineutral). (c,d) Ion energy flux 
spectra [in D units] of pure plasmas (solid), with carbon (dashed), with deuterium 
(dash-dotted), and with carbon and deuterium impurities, (c): main ion hydrogen, 
(d): helium, (e) and (f): carbon and deuterium particle flux spectra. 



impurities (dotted); ric/rii = 0.04167 corresponding to Z cff = 2 in hydrogenic plasmas, 
and ri£)/ni = 1/3. The stabilizing effect of carbon on ITG modes, seen also in the linear 
simulations, appears as a strong reduction (by 51%) of the ion energy fluxes that acts 
at all wave numbers (we note that TE modes are stabilized by collisions). The presence 
of deuterium slightly reduces the transport at higher wave numbers (this feature is 
consistent with corresponding linear simulation results, not shown here). Interestingly, 
when there is carbon impurity in the plasma, the presence of a deuterium minority does 
not reduce the absolute value of the transport, only the spectrum gets shifted slightly 
towards lower wave numbers. The effect of these impurities is qualitatively the same 
on the electron energy flux, the difference is that the reduction in the electron energy 
transport is not as strong as in the ion energy transport (only 35%). Figure [TJ (d) 
shows how the presence of D and C affects the energy transport in a helium plasma. 
The effect of a C impurity is considerably weaker in a helium plasma than in hydrogen 
isotope plasmas and it is concentrated to lower wave numbers. The fact that we kept 
ric/ni fixed in the simulations [shown in Figs. [T] (c) and (d)] instead of n c /n e is part 
of the reason why we find a weaker effect of carbon in the simulation with helium than 
with hydrogen main ions, since then carbon density is smaller in the helium plasma. 
However, even if we kept nc/n e fixed, looking at Fig. [T] (a) we expect to have a weaker 
impact of the carbon impurity on transport because the relative change in the linear 
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growth rates (compared to the pure plasma case) for given carbon concentration is 
smaller in a helium plasma than in hydrogen. While both carbon and deuterium act to 
reduce the energy transport in the helium plasma (compare the gyro-Bohm normalized 
fluxes in the legend of Fig. Q] (d)), interestingly, when both impurities are present the 
absolute value of the transport remains almost the same as in the pure plasma, however 
the shape of the energy flux spectra changes somewhat, being more peaked around its 
maximum. 

In Figs. Q] (e,f) the carbon and deuterium particle fluxes are shown in impure 
hydrogen and helium plasmas, respectively (positive sign corresponding to an outward 
particle flux). In species units the magnitude of the carbon particle flux is the same in 
a hydrogen and a helium plasma when the ratio nc/rii is kept fixed; compare the solid 
and dash dotted curves in Fig. [T] (e). The shape of the carbon particle flux spectra are 
similar to the ion energy flux spectra in the hydrogen and helium plasmas [H+C and 
He+C curves in Figs. (c,d) compared to those in Fig. (e)], however, the magnitude 
of Qi is smaller in the H+C than in the He+C plasma. It is also interesting to note, 
that in the hydrogen plasma the carbon particle flux remains almost unaffected when a 
deuterium minority is introduced; in contrast to the helium plasma, where the carbon 
flux drops dramatically (by 63%) in the presence of deuterium; a trend opposite to what 
was found for ion energy fluxes. 

Figure Q] (f ) shows that the deuterium particle fluxes in species units appear to be 
almost three times higher in a helium plasma than in a hydrogen plasma, even if the 
nr>/rii values are held constant, so the number of deuterium particles is half as many 
in the helium plasma as in the hydrogen plasma. In spite of the fact that the linear 
growth rates and the energy transport are reduced when carbon impurity is present, the 
deuterium particle fluxes increase both in a H(+D) and a He(+D) plasma when carbon 
impurity is introduced. 

3. Impurity peaking factor in the presence of poloidal asymmetries 

In this section we calculate the zero flux density gradient (peaking factor) of trace 
impurities under the effect of a poloidally varying equilibrium electrostatic potential 4>e 
that can arise due to e.g. the presence of a species with strong temperature anisotropy 
[20] . Especially, in case of RF heating on the outboard side, the trapped population of 
the heated particle species increases in the vicinity of the resonance position [2Tj ITB] . 
The peaking factor depends on the linear mode characteristics and on the form of <pE- 
Here, we neglect effects caused by a radial electric field, such as toroidal rotation. 

We assume eA<pE/Tj <C 1, where A<pE refers to the poloidal variation of 4>e, so 
that poloidal asymmetry effects on the main species can be neglected; this assumption 
is needed to justify the use of GYRO simulations to obtain the linear mode characteristics. 
However, we allow for ZeA(f>E/T z = 0(1). 
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3.1. Perturbed impurity distribution 

We consider the particle transport driven by a single, representative, toroidal mode. 
The impurity peaking factor is calculated by requiring the linear impurity flux T z = 
S {—k y h z (f)* I B) = Q (—k y J d 3 vJQ(z z )g z (f>*/B^ to vanish (steady state is assumed and 
impurity sources are neglected). Here (•) denotes a flux surface average, h z is the 
perturbed impurity density, g z the non-adiabatic part of the perturbed impurity 
distribution function, J is the Bessel function of the first kind, z z = k±v±/u cz , 
k± = (1 + s 2q d 2 ) 1 ^ 2 k y , and the rest of the notation is standard. The species indexes 
z, e and i refer to impurities, electrons and main ions. 

The linearized GK equation is to be solved to obtain the non-adiabatic part of the 
impurity distribution function g z , 

- i{u - u Dz - u E )g z - C[g z ) = -i^j^- (u - ujJ z ) <f)J (z z ), (1) 

where $ is the extended poloidal angle, u = u r + 27 is the mode frequency, f z o(i>, £) = 
n z o(?p){m z /2iiT z (tp)} 3 / 2 exp{— E /T z (ip)] is the equilibrium Maxwellian distribution. \i = 
m z v 2 J (2B), 8 = m z v 2 /2 + Ze<\> E , a£ = u* z [1 - L nz Zed r {^) E e(j) E /T z - 3/2)L nz /L Tz ] 
with u* z = -k y T z /ZeBL nz , L nz = - [d r (lnn z )] _1 and L Tz = - [d r (hiT z )]~ . 
Furthermore, ud z = —2k y /m z (£ — Ze<p E — fiB /2)V('d) /u cz R, where Vi^d) = costf + 
si)sin'd. The E x B drift frequency in the equilibrium electrostatic field is u E = 
k y [d r (j) E — s$d$((j) E )/r} /B. Our ordering Ze<f> E /T z ~ 0{\) requires that u E /u is 
formally ~ 1/Z small; this requirement is mostly satisfied for experimentally relevant 
poloidal asymmetries. We consider ion scale modes, Z{ < 1, and we allow the parallel 
streaming term of main ions to be comparable to the mode frequency. In the following we 
perform a perturbative solution of the impurity GK equation ([1]) in the small parameter 
Z~ x l 2 <c 1. We assume that impurity self-collisions dominate over unlike-particle 
collisions, requiring n z Z 2 /n e = 0(1), thus we model only impurity- impurity collisions 
with the full linearized collision operator C z } [•] . 

Keeping in mind that u E /u, ud z /u, ojJ z /ou, and Jq{z z ) — 1 ~ — z 2 /4 are all 
~ 1/Z small, the lowest order equation —iugo — C z }[go] = —iuZe<j)f z o/T z for g z = 
9o + 9i + 92 + • • • is satisfied by g = Ze<f)f z0 /T z , since C zz [g oc f z0 ] vanishes. To first 
order in Z~ 1/2 , we have v\\d#(g Q ) / '(qR) - iug x - C { z ] z [gi] = 0. Using Cfj^i oc v\\f z0 ] = 0, 
and de{fzo)\e = 0, we obtain g\ = —iZef z oV\\d$(4))/(T z uqR). To next order, the GK 
equation reads 

- iug 2 - C ( p z [g 2 ] = -i(uj Dz + uj e ) 9o - v\\d#( gi ) / (qR) + iZe<j)f z0 (ujz 2 /4 + ujJ z )/T z . (2) 

Using that m z v\\d^(v\\)\s^ = —fid^B — Zed$(j) E , the parallel compressibility term in the 
right hand side of Eq. §2§ , can be expressed as 

v^dgi = . Zef z0 f 2 d_ / }_df\ _ }_d(p 
qR d'd ~ l uT z q 2 R P 1 d$ \Rd$) Rd$ 



n dg z 
qR dd 



vj dlnB Ze d<p E 
2 d'd m z d'd 



(3) 
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We decompose the magnetic drift frequency and the FLR parameter as ud z = ojdx+^d\\ 
and z 2 = z 2 —z 2 , where ojd x and z 2 x are proportional to v 2 and lod\\ and z 2 are proportional 
to v 2 . Also, we write g% = §2 + §2, where §2 has contributions proportional to f Z Q and 
v 2 f z0 , thus C z }[g 2 ] = 0. The collisionless, isotropic part of Eq. (J2J) is solved by 



92 = go 



Zecj) uz 2 JA + q& 

■JzO r 



-.0 



co 2 T z q 2 R 2 d$ 



v 2 d\nB Ze 



2 <9$ 



(4) 



Then (72 should satisfy the remaining part of Eq. (j2J) which, upon division by —iu reads 

v 2 7„„,2j 



#2 



Co' 



#0 



Ze0 ^ 



II 





"1 df 







+ 



<91nB 



2R&& d$ 



(5) 



Equation [5] can be written as g 2 — [fife] ~~ -^ u y/zo = 0, where X is independent of 
velocity. The solution can be written in the form g 2 = H(v)/ z q, where 5 can incorporate 
any nontrivial velocity dependence. Taking the density moment of Eq. (J5]) 



d 3 v[g 2 -iCU[g 2 ]/uJ-Xvlf zQ ] = 
= J d 3 vE(v)C[ l) z [f z0 ] = 0, Eq. 



0. (6) 
(HO) can be reduced to 



Since Jd 3 vC ( J>[E(v)f z0 ] 
fd 3 v[g2-Xv 2 f z0 \ =0. 

Since we are only interested in the density moment of the perturbed distribution 
function needed to evaluate the impurity peaking factor, instead of solving the 
complicated collisional equation <^ unnecessarily, we derive a substitute function g 2s 
that has the same density moment as g 2 ; we choose g 2s = Xv 2 f z0 /3. 

We introduce a substitute function for g z , denoted by g s , with the property 
/ d 3 vJ (z z )g s = J d 3 vJ (z z )g z + 0(Z' 2 )\ we define it as g s = g + g 2 + g 2s , that reads 



9 S T Z 
Zef, 



-.0 



4u Dx /3 + uj e - uj z _ z; 
uj 6 



21 



3(uqR¥ 



x. 



,d 2 (j> 3d<pd_ ( Ze<p E 



(7) 

where x z = v/v z , with v z = (2T z /m z ) l l 2 . It is important to emphasize, that g s is not 
an approximate solution to Eq. (JTJ), but merely a function that has the same density 
moment as the solution. Taking other velocity moments of this function would give 
erroneous results, and g s cannot be used to proceed with the perturbative solution. 



3.2. Zero flux impurity density gradient 

To calculate the impurity peaking factor we assume a simple sinusoidal poloidal 
asymmetry of the form Ze<pE/T z = —kcos(6 — 5) where k is the asymmetry strength 
and S determines the poloidal position of the impurity accumulation. When evaluating 
the density moment of Eq. (j7|), we keep the 0{Z~ 1 ) correction from Jq{z z ) where it 
multiplies go- Also, we neglect 0(e) corrections when solving for a/L^ z , the value of 
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a/L nz where (T z ) = 0, to find 



2aVi Zrrii u r 
(qR ) 2 k y pi m z lo 2 + i 1 



06 




where pt = Vimi/eB , and (...)<£ = (. . . J\f\(j)\ 2 ) / (J\f\(j)\ 2 ) , with A/" = exp[«;cos(0 — 6)]. 
The FLR terms do not appear in a/L° nz , since their imaginary part is zero. The 
first and second terms of Eq. (jHJ) represents the contributions from ujd and wg, 
respectively. Since the last term of Eq. (jHJ) contains only non-negative quantities, 
except u r , impurity parallel dynamics acts to increase(decrease) the impurity peaking 
if bj r is negative(positive). Note, that we use the sign convention of GYRO; u r is 
negative for modes propagating in the ion diamagnetic direction. The impurity parallel 
compressibility term introduces a dependence on mode frequency, and also on the charge 
to mass ratio, consistently with Ref. 




Figure 2. (a) Shear dependence of the impurity peaking factor. Thin solid line 
represents the poloidally symmetric case; the corresponding gyro results are shown 
with red diamonds. Thick solid line represents an in-out asymmetry, k = 0.5. The 
contributions from the different terms in Eq. ((5J) arc also plotted; dashed line: magnetic 
drifts, dash-dotted: E x B drifts, dotted line: parallel compressibility, (b) Shear 
dependence of uj r and 7 from GYRO. 



The shear dependence of the peaking factor of a fully ionized Nickel impurity 
calculated from Eq. (JHJ) is shown in Figure |2Ja). The local profile and geometry 
parameters assuming a circular model Grad-Shafranov equilibrium are: r/a = 0.3, 
R /a = 3, k e p s = 0.3, q = 1.7, a/L ne = 1.5, Tj/T e = 0.85, a/L Te = 2 and a/L Ti = 2.5, 
and collisions are switched off in the simulations to avoid discrepancies due to the 
different collision operators in this paper and in gyro. The poloidally symmetric case 
is shown with thin solid line and shows a relatively good agreement with the exact values 
obtained from gyro simulations represented by the diamond markers; the differences 
towards higher magnetic shear appear as a result of an overestimation of the parallel 
compressibility term due to the approximations done in the derivation of Eq. (jHJ). A 
scenario with in-out impurity asymmetry for the asymmetry strength k = 0.5 is also 
plotted with thick solid line, together with the contributions from the different terms 
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in Eq. (JHJ). Around s = the contribution from ue (dash-dotted line) vanishes, and in 
this specific case the parallel compressibility term (dotted line) is also small due to the 
small value of u r ; the peaking factor is dominated by u>d (dashed line). If |w r |/7 <C 1, 
that happens for lower shear in the case studied here [see Fig. [2] (b)], the frequency 
dependent part of the third term in Eq. (jHJ) is small, thus parallel compressibility effects 
do not play an important role; see the dashed curve in Fig. [2] (a). 

As the magnetic shear increases the ue term starts to reduce the peaking factor, 
since (9 sin(6 l — n))^ is negative for moderately ballooning eigenfunctions. In Eq. (jHJ) the 
shear dependent part of T>, and the ue term can be combined as (as / Rq) (9 sin 9) ^(2 ± 
ft/e), where e = r/Ro, and the plus(minus) sign applies for out-in(in-out) impurity 
asymmetry. Thus, if (f)(9) and u are only weakly dependent on shear, increasing shear 
is expected to reduce the impurity peaking in the in-out asymmetric case when 2 < n/e 
[for our simulation parameters n/e = 5]. From Eqs. (10) and (12) in Ref. [19] one can 
estimate the ratio of the asymmetry strength and the inverse aspect ratio; assuming 
T e ~ Tj ~ T 2 , and taking the limits ^> 1 and e < 1 so that eax is also small we 
obtain n/e « a T ZA" m /(l + Z cS ), where ay is the minority temperature anisotropy, and 
X m is the minority fraction. 

4. Conclusions 

In the present paper we first study how impurity transport is affected by the charge and 
mass of the main ion species and the ion composition, and what impact the impurities 
have on energy transport, then we consider the effect of a poloidally varying equilibrium 
potential on impurity peaking. 

Linear gyrokinetic analysis shows that for increasing impurity concentration the 
growth rates change in a linear fashion; in main ion species units the change in the 
growth rates when increasing impurity concentration to a given value - i.e. the slope of 
the ^ya/c S i(nc/n e ) curves - is almost independent of the charge and mass of the main 
species. In the studied case, based on local profile and geometry data in a deuterium 
discharge from DIII-D, for ITG modes carbon impurity has a strongly stabilizing,- for 
TE modes (when neglecting collisions) a weak destabilizing - effect. For TE modes 
the effect is almost purely due to plasma dilution, for ITG modes the presence of the 
impurity itself has a comparable effect to the dilution that still dominates. The strong 
stabilizing effect of carbon in hydrogenic plasmas appears in nonlinear simulations as a 
strong reduction of energy transport, the relative change in the energy fluxes being even 
higher than that of the linear growth rates. The effect of impurities on energy transport 
is smaller in a helium plasma, due to the higher species units growth rates and the 
lower main ion concentration. The presence of deuterium mainly just changes the shape 
of the energy flux spectrum, although slightly reduces the transport in a hydrogen 
plasma. In general the presence of impurities, even in relatively small quantities, can 
cause significant deviations from a pure gyro-Bohm mass and charge scaling of energy 
fluxes. (Again, by "pure gyro-Bohm scaling" we refer to the situation when differences 
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in the turbulence between similar plasmas with different ion composition can be exactly 
transformed out by normalizing temporal and spatial scales to a/vi and pi, respectively. 
We do not mean that the transport is not gyro-Bohm in the general sense.) 

We conclude that, while the apparent deviations from a pure gyro-Bohm scaling, 
- due to charge, electron-to-ion mass ratio or collisional effects, and importantly due 
to the presence of impurities, - can mostly be explained and understood from a linear 
analysis of the underlying microinstabilities, impurity transport do not follow these 
naive expectations. In particular, if we would normalize the impurity fluxes to ion 
energy fluxes we would get significantly different values in plasmas with different main 
ion species. Thus we might expect that the steady state impurity profiles would also vary 
in these plasmas (especially when impurity sources and neoclassical impurity transport 
are accounted for), causing further deviations from a pure gyro-Bohm scaling. 

A poloidally varying electrostatic field, appearing due to e.g. a poloidally 
asymmetrically distributed RF heated particle species, even if being too weak to modify 
the dynamics of the main species, can essentially change impurity transport compared 
to a poloidally symmetric situation. We show that the combined effect of the arising 
poloidally asymmetric impurity distribution and the E# x B^, drift of impurities can 
lead to a reduction or even a sign change in the impurity peaking factor. The effect 
becomes important at high impurity charge, when the magnetic and diamagnetic drifts 
(proportional to 1/Z) become as small as the E# x B^ drift in the poloidally varying 
equilibrium potential. We demonstrate that to lowest order in 1/Z finite Larmor 
radius effects and impurity collisions do not affect impurity peaking driven by ion scale 
microinstabilities as long as the impurity collisions are dominated by self-collisions and 
the impurity collision frequency is not much larger than the mode frequency. We present 
and analyze a simple analytical expression for the impurity peaking factor including 
contributions from the Eg x B^, drift and parallel compressibility, and depending on the 
linear mode characteristics. We find that to get sign change in the impurity peaking 
factor, a necessary criterion is that the ratio of the asymmetry strength k and the inverse 
aspect ratio e is larger than two. 
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